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Abstract 

We derive the distribution of the center of mass S of the integrated 
superBrownian excursion (ISE) from the asymptotic distribution of the 
Wiener index for simple trees. Equivalently, this is the distribution of the 
integral of a Brownian snake. A recursion formula for the moments and 
asymptotics for moments and tail probabilities are derived. 
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1 Introduction 

The ISE (integrated superBrownian excursion) is a random probability measure 
on E"*. The ISE was introduced by David Aldous P as an universal limit object 
for random distributions of mass in M.'^: for instance, Derbez & Slade JT] proved 
that the ISE is the limit of lattice trees for d > 8. The ISE can be seen as 
the limit of a suitably renormalized spatial branching process (cf. |HE1)) or 
equivalently, as an embedding of the continuum random tree (CRT) in W'-. The 
ISE is surveyed in [T3 . 

Formally, the ISE is a random variable, denoted J^, with value in the set of 
probability measures on R'^. In Section 13 we give a precise description of J in 
terms of the Brownian snake, following |14l Ch. IV. 6]. As noted in even in 
the case d = 1, where the support of J is almost surely a (random) bounded 
interval denoted [R, L] , little was known about the distributional properties of 
elementary statistics of such as R, L, or the center of mass 



J xj{dx) 
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This is still true, but recently R — L was shown |S] to be the law of the radius 
of some model of random geometries called fluid lattices in Quantum geometry 
or random quadrangulations in combinatorics, and the law of R and L were 
investigated by J.F. Delmas 8 . The joint law of {R,L) is still unknown, as far 
as we know. 

The Wiener index w{G) of a connected graph G with set of vertices V{G) is 
defined as 

wiG) = ^ d{u,v), 

{u,v}<£V{G) 

in which d{u, v) denotes the distance between u and v in the graph. For G a 
random simple tree with n nodes, w{G), suitably normalized, is asymptotically 
distributed as ^—r], in which ^ and i] are the following statistics of the normalized 
Brownian excursion (e(f))o<t<i: 

^ = 2 [ e{t) dt (1) 

77 = 4 / min e{u) ds dt (2) 

Jo<s<t<i 

(cf. |12j. where the joint moments of (77,0 computed from explicit formulas 
holding for random binary trees). 

In this paper, we derive asymptotics, and explicit induction formulas, for 
the moments of S, from similar results for 77. 



2 The ISE and the Brownian snake 

We recall briefly the description of J' in terms of the Brownian snake, from 
[H Ch. IV.6] (see also [TOl HHI)- The UfeUme process ( = (C(O)teT of 
Brownian snake is a stochastic process with values in [0, +00). Let the Brownian 
snake with lifetime C be denoted 

W = (W-.(i))o<.<i, o<t<c(.) ■ 

The Brownian snake can be seen as a description of a "continuous" population 
T, through its genealogical tree and the positions of its members. The lifetime C 
specifically describes the genealogical tree, and W describes the spatial motions 
of the members of population T. 

A member of the population is encoded by the time t it is visited by the 
contour traversal of the genealogical tree, ({t) being the height of member t ^ T 
in the genealogical tree (C(i) can be seen as the "generation" t belongs to, or 
the time when t is living). Let 

C(s,i)= min C(")i s A i = argmin C(w). 

s<u<f s<u<t 
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Due to the properties of the contour traversal of a tree, any element of s A t is a 
label for the more recent ancestor common to s and t, and the distance between 
s and t in the genealogical tree is 

dis,t) = as)+at)-2C{s,t). 

If it is not a leaf of the tree, a member of the population is visited several times 
(fc + 1 times if it has k sons), so it has several labels: s and t are two labels of 
the same member of the population if d{s, t) = 0, or cquivalcntly if s At D {s,t}. 
Finally, s is an ancestor of t iff s € s A t. In this interpretation, Ws{u) is the 
position of the ancestor of s living at time u, and 

Ws = Ws (C(s)) , s e T, 

is the position of s. Before time m = C{si,S2), si and S2 share the same 
ancestor, entailing that 

iWsAt))o<t<m = iWsAt))o<t<m- (3) 

Obviously there is some redundancy in this description: it turns out that the 
full Brownian snake can be recovered from the pair (Ws,(^{s)] (see QHl 

V /0<s<l 

for a complete discussion of this). 

In the general setting |14[ Ch. IV], the spatial motion of a member of the 
population is any Markov proces with cadlag paths. In the special case of the 
ISE, this spatial motion is a d-dimensional Brownian motion: 

a) for all < s < 1, t ^ Ws{t) is a standard linear Brownian motion started 
at 0, defined for < t < C(s) ; 

b) conditionally, given C, the application s W^s(-) is a path- valued Markov 
process with transition function defined as follows: for si < S2, condi- 
tionally given Wsj(.), (W^«2 ("^ + ^))o<t<c(s2)-m ^ standard Brownian 
motion starting from Wsj(m), independent of Wsi(-)- 

The lifetime C is usually a reflected linear Brownian motion |14) , defined on 
T = [0, -|-oo). However, in the case of the ISE, 

C = 2e, 

in which e denotes the normalized Brownian excursion, or 3-dimensional Brow- 
nian bridge, defined on T = [0, 1]. With this choice of the genealogical tree is 
the CRT (see ^), and the Brownian snake can be seen as an embedding of the 
CRT in M.'^. We can now give the definition of the ISE in terms of the Brownian 
snake with lifetime 2e 14, Ch. IV. 6]: 

Definition 2.1. The ISE J' is the occupation measure of W. 



3 



Recall that the occupation measure J of a. process W is defined by the 
relation: ^ 

f{x)J{dx) ^ j^f {W,) ds, (4) 
holding for any measurable test function /. 

Remark 2.2. It might seem more natural to consider the Brownian snake with 
lifetime e instead of 2e, but we follow the normalization in Aldous pp. If we 
used the Brownian snake with lifetime e instead, W , J and S would be scaled 
by 1/^2. 

3 The basic identity 

Theorem 3.1. Let N be a standard Gaussian random variable, independent of 
rj. Then 

S '= ^ N. 

Proof. Based on the short account of the Brownian snake theory in Sectional the 
proof is now pretty easy. Specializing I^J to f{x) = x, we obtain a representation 
of 5: 

[ Ws ds, 







which is the starting point of our proof. We have also, directly from the defini- 
tion of the Brownian snake. 

Proposition 3.2. Conditionally, given e, [Ws] is a Gaussian process 

\ / Q<s<l 

whose covariance is C{s,t) = 2mins<t(<t e{u), s <t. 

Proof. With the notation ( = 2e and m = C{si,S2) — 2miiis-^<cu<s2 we 
have, conditionally, given e, for si < S2, 

Cov(Ws„Ws,) = Cov{W,,{asi)),WsMs2))-Ws,{m) + WsAm)) 

= CoY{W,,{asi)),WsAm)) 
= Cov(W^,,(C(si)),M^.,(m)) 
= m, 

in which b) yields the second equality, Q yields the third one, and a) yields the 
fourth equality. □ 

As a consequence of Proposition 13.21 conditionally given e, 5 is centered 
Gaussian with variance 



C(s, t) ds dt = rj. 

'[0,1V 

This last statement is equivalent to Theorem 13. II □ 

Remark 3.3. The d-dimensional analog of Theorem 13.11 is also true, by the 
same proof. 
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4 The moments 



We know that the odd moments of S vanish. 
Theorem 4.1. For k>0, 

^L* J - 2(9^-4)/2r((5fc-l)/2) 
in which is defined by ai ~ I, and, for k > 2, 

fe-i 

flfc = 2(5fc - 4)(5fc - 6)afe-i + ^ a^ak-i- (6) 

Proof. From Theorem 13.11 we derive at once that 

E [S^''] = E [r]''] E [AT^fe] . 
As a special case of ^1 Theorem 3.3] (where ak is denoted ujQf.), 

^ = 2(^^-4)/2r((I-i)/2) 

and the result follows since E [iV^fcj ^ (2fc)!/(2'=fc!). □ 
In particular, see again 12, Theorem 3.3], 

E [S"^] = E [?/] = 

E [5"*] = 3E [rf] = 7/5, 

as computed by Aldous 1, by a different method. (Aldous' method extends to 
higher moments too, but the calculations quickly become complicated.) 
We have the following asymptotics for the moments. 

Theorem 4.2. For some constant f3 = 0.981038 ... we have, as k — > oo, 

afc~/3 50^-1 (fc-l)!2, (8) 

E [r;'^] ^ ^f^k'/^ (5e)-'=/2 k'^/^ (9) 



E [S^"] ^ ^^!^!^(2fc)i/2 (lOe^)"''/' (2fc)i•2^ (10) 

Proof Set 



50'^-i(/c-l)!2' 

We have 6i = ai = 1 and, from ©, 62 = ^3 = 49/50 = 0.98. For fc > 3, © 
translates to 



6. = 1 - 



1 \ , J2i=l f^if^k-i 



25(fc-l)V ^-'"^ ' 50'=-i(fc-l) 



Ek-2 
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Thus bk increases for fc > 3. We will show that bk < I for all fc > 1; thus 
(3 = limfe^oo bk exists and JSjl follows. Stirling's formula and 10, Q then yield 
and (O. 

More precisely, we show by induction that for fc > 3, 

bk<l- ^ 



25(fc- 1)' 



This holds for = 3 and fc = 4. For fc > 5 we have by the induction assumption 
bj < 1 for 1 < j < fc and thus 



Sk 



Ek — 2 
i—2 ^i^k — i 

50*^-1(^-1)!^ 



(11) 



_ 1 ^^ (z^l)P(fc-z-l)P ^^ 
"50^ (fc-l)!2 ^'""-^ 

< \ 2 (2 + (fc-5) ^1 

" 50(fc-l)'(fc-2)' V [k-SfJ 

1 

< 9 • 

25(fc-l)(fc-2)^(fc-3) 

The induction follows. Moreover, it follows easily from Hll() that bk < P < 
6fc + ^(fc-2)-3. 

To obtain the numerical value of /3, we write, somewhat more sharply, where 
< 61 < 1, 

1 4 

= ■ 79— -n&2&/c-2 H ■ 3— 79- —bzbk-a 

25 (fc - 1)^ (fc - 2f 25 (fc ~ 1)^ (fc - 2f (fc - 3)2 

3fi 

+ ^(k - 7) 2 2 

25 (fc - 1)^ (fc - 2)^ (fc - 3)2(fc - 4)2 

and sum over k > n for n = 10, say, using < bk^2 < P and 5„-2 < ^fe-a < P 
for k > n. It follows (with Maple) by this and exact computation of &i, . . . , biQ 
that 0.981038 < /3 < 0.9810385; we omit the details. □ 

Remark 4.3. For comparison, we give the corresponding result for ^ defined 
in 0. There is a simple relation, discovered by Spencer and Aldous [2], 
between its moments and Wright's constants in the enumeration of connected 
graphs with n vertices and n + k edges .19; , and the well-known asymptotics of 
the latter lead, see ^1 Theorem 3.3 and (3.8)], to 

E[C''] ~3V2fc(3e)-'=/2fcfc/2. (12) 
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5 Moment generating functions and tail esti- 
mates 



The moment asymptotics yield asymptotics for the moment generating func- 
tion E [e*'^] and the tail probabilities ¥ {S > t) as t ^ oo. For completeness 
and comparison, we also include corresponding results for rj. We begin with a 
standard estimate. 



Lemma 5.1. (i) If ^ > and & e R, then, as x oo 

oc 
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a; ^ I — I [e 'xj e 

fc=i 



(ii) // —oo < 7 < 1 and 6 G K, then, as x — > oo, 

^''^'^^^^ ~ (1 - ^)-i/2(g7a,^''/(i-T)e(i-7)(e-'x)i''<i-^' 

fc=i 

The sums over even k only are asymptotic to half the full sums. 

Sketch of proof . (i). This is standard, but since we have not found a precise 
reference, we sketch the argument. Write k^k~''^x^ — e^^^^ — fc^'e^'-'^^ where 
g{y) — —72/ In 2/ + ?/ In a; and /(y) — g{y) + b\n.y. The function g is concave with 
a maximum at yo = 2/0(2;) = e^^x^/'^ . A Taylor expansion yields 



k=l 







1/2 
Va 



e-V2,, 

V 7 / 



(ii). Follows by (i) and Stirling's formula. □ 

Combining Theorem 14. 21 and Lemma l5. II we find the following asymptotics 
for E [e*"] = Efc E [r]^] t^/k\ and E [e*'^] ^ Y.k cven^ [S^] t^/k\. 



Theorem 5.2. As t 00, 

E[e"]~^*.""°, (13) 

9l/2_3/2o 

E [e*^] ~ J/, ^ t^e*V40. (14) 

Proof. For 77 we take & = 1/2, 7 = 1/2, x = (Se)^^''^/; in Lemma fS.lf iil: for S 

we take 6 = 1/2, 7 = 3/4, x = (10e3)-i/4t_ □ 
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The standard argument with Markov's inequality yields upper bounds for 
the tail probabilities from Theorem 14 . 21 or 

Theorem 5.3. For some constants Ki and K2 and all x > \, say, 

P(?7 > < ii:ia;exp(-|a;2), (15) 
I'{S>x)< i^2a;2/3exp(-|lOi/^a;''/^). (16) 

Proof. For any even fc and a; > 0, P(|5| > x) < x'^W. [S''=] . We use (HHJ and 
optimize the resulting exponent by choosing k = lO^/^x'*/^, rounded to an even 
integer. This yields 1)16(1 : we omit the details. 115(1 is obtained similarly from 
(jnj, using k = [5x^J . □ 

Remark 5.4. The proof of Theorem 15.31 shows that any Ki > 27r^/^/3/5^/-^ w 
4.9 and K2 > 10^/^/37r^/^/5 w 1.6 will do for large x. Alternatively, we could 
use Theorem 15.21 and ¥ {S > x) < e~*'^E [e*'^] for t > and so on; this would 
yield another proof of Theorem 15.31 with somewhat inferior values of Ki and 
K2. 

The bounds obtained in Theorem 15.31 are sharp up to factors 1 + o(l) in 
the exponent, as is usually the case for estimates derived by this method. For 
convenience we state a general theorem, reformulating results by Davies [2| and 
Kasahara [Oj . 

Theorem 5.5. Let X be a random variable, let p > and let a and b be positive 
real numbers related by a ^ \/(j)eW) or, equivalently, b — {pea)^^/P . 

(i) IfX>0 a.s., then 

- \nF{X > x) r^axP as x 00 (17) 

is equivalent to 

(EX^y/'' - fer^/P as r-^ 00. (18) 

Here r runs through all positive reals; equivalently, we can restrict r in 
((18(1 to integers or even integers. 

(ii) If X is a symmetric random variable, then ((17() and p8(l are equivalent, 
where r in ((18(1 runs through even integers. 

(iii) If p > 1, then, for any X , ((17() is equivalent to 

ln(Ee*^)-ci« as t ^ 00, (19) 

where l/p+ l/q = 1 and c = q~^{pa)^^'^^^^ — q^^e'^^^W . (This can also 
be written in the symmetric form (paY^qcY = 1 and as b — e^^^^'{qcY/'^ .) 
Hence, if X > a.s., or if X is symmetric andr restricted to even integers, 
((19(1 is also equivalent to ((1811 . 
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Proof. For X > 0, (i) and (iii) are immediate special cases of more general 
results by Kasahara |13[ Theorem 4 and Theorem 2, Corollary 1], see also |5| 
Theorem 4.12.7]; the difficult implications ifTH)) =^ ((T7jl and =» ((T7|l 
were earlier proved by Davies [7j (for p > 1, which implies the general case of 
(i) by considering a power of X). Moreover, (|19|l =^ (|17|) follows also from 
the Gartner-Ellis theorem P| Theorem 2.3.6] applied to n~^^PX. Note that, 
assuming X > 0, (EX^Y^^ is increasing in r > 0, which implies that H18() for 
(even) integers is equivalent to H18I) for all real r. 

(ii) follows from (i) applied to \X\, and (iii) for general X follows by consid- 
ering max(X, 0). □ 

We thus obtain from Theorem 14. 21 or I5T^ the following estimates, less precise 
than Theorem 15 . 31 but including both upper and lower bounds. 

Theorem 5.6. As x cx3, 

ln{¥iri>x)) ^-^x^, (20) 

ln(P (S > x)) - - JlO^/^ x^^\ (21) 

Proof. For rj we use Theorem 15.51 with p = q = 2, b = {5e)~^^'^, a — 5/2 and 
c = 1/10. For S we have p = 4/3, q ^ A, b = (lOe^)-'^/*, a = 3 • 10^/^4 and 
c=l/40. □ 

Remark 5.7. I|20|l can also be proved using the representation (0) and large 
deviation theory for Brownian excursions, cf. §5.2] and ^Dl- The details may 
perhaps appear elsewhere. 

Remark 5.8. (|21() can be compared to the tail estimates in ^ for the density 

function of Wu, the value of the Brownian snake at a random point U, which 
in particular gives 



In 



{¥(Wu > x)) 3-2-5/3 2;4/3_ 



Remark 5.9. For ^ we can by (|12|l use Theorem 15.51 with p = q = 2, b ^ 
(3e)-i/2, a = 3/2 and c = 1/6. In [12, the variable of main interest is neither ^ 

nor 7? but C C-V- By Minkowski's inequality E [^'=] < E [(*=] ^^''+E [rj''] ^'^ 
and ®, 113 follows 

1 1 (Efc'^])'^' mc^X)^"' 1 

< lim inf ^ < lim sup < 



/3i V5i fe^oo ^.^^ /g^- 



This leads to asymptotic upper and lower bounds for ln(P (C > a;)) too by [7] or 

ini. We can show that limfe^oo A:"^/^ (E [C''])^^'' and lim^^^oo a;~Mn(P (C > x)) 
exist, but do not know their value. 
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6 Concluding remarks 



The center of mass of the ISE turns out to be related to the Wiener index 
of simple trees, but note that S is also related to the Wiener index of ran- 
dom planar quadrangulations: let (&n,e„)) denote the uniform choice of a 
quadrangulation with n faces, and of a "marked" oriented edge in it, and set 

As a consequence converges weakly to c • {S ~ L), where L is 

the left endpoint of the support of J^, and c is a known constant. The joint law 
of (5*, L) is not known, as far as we know. 

Acknowledgement. We thank Jim Fill for helpful comments. 
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